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VELOCITY DISTRIBUTION IN THE BOUNDARY LAYER 
OF A SUBMERGED PLATE.* 
By 14. Hansen, 

This report deals with the measurement of the velocity dis- 
tribution of the air in the vicinity of a plate placed parallel 
to the air flow. 

I. Notation, Definitions and Theoretical Results 

p, density of the air, 

p., viscosity of the air, 

v, .kinematic viscosity coefficient, 

x, distance of test point from leading edge of plate, 

y, distance of test point from plate, 

I, length of plate, 50 cm (19.69 in.), 

b, width of plate, 38 cm (14.96 in.), 

U, air velocity outside of boundary layer, 

u, component of boundary layer air velocity paxallel to 
plate, 

5, thickness of boundary layer, 

t 0 , shearing stress on plate, 

♦"Die Geschwindigkeitsverteilung in der Grenzschicht an einer 
eingetauchten Platte." From Abhandlungen aus dem Aerodynami- 
schen Institut an der Technischen Kochschule Aachen, No. 8^ 1928, 
pp. 31-45. 
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R T = L characteristic coefficients. 

L. Prandtl and H. Blasius succeeded in determining theoret- 
ically the velocity distribution for the case of a lajninar flow.* 
Their results may be briefly summarized as follows. 

Let x and y be the coordinates, respectively, parallel 
and perpendicular to the plate, and U the velocity of the un- 
disturbed air flow. The velocity curves are similar for differ- 
ent values of x, the distance at which the same velocity oc- 
curs, increasing with the square root of x. The velocity 
obeys the law: 

For comparison with experimental results, it is convenient 
to introduce a "boundary-layer thickness," although no exact 
value can be assigned to it, because the velocity u changes 
asymptotically to U. For example, y may represent the boun- 
dary-layer thickness for which the function f differs from 
unity by a certain fraction (say If the function f in 

equation (l) is approximated by a parabola, which changes asymp- 
totically into the line u = U, then 

*H. Blasius, Zeitschrift fur Mathematik und Physik, Vol* 56, No, 
1, 1908. 
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is then obtained for the transition point. This value of y 

will be taken for comparison with the experimental results. 

The surface friction is represented by t q = u (4—^) 

• y y=o 

According to the Blasius theory, it is 

t = 0.332 p U 2 — (3) 

.Various observations led Prandtl to surmise that the lami- 
nar flow changes, with increasing Reynolds Number, to a turbu- 
lent flow, similar to the flow through a tube. Since the ele- 
mentary law of turbulent friction is not yet known, the veloc- 
ity distribution cannot be theoretically determined for this 
case. Nevertheless, if the velocity distribution experimentally 
determined for tubes and the likewise experimentally determined 
law of resistance be transferred to the case of the boundary 
layer-, the thickness of the boundary Layer can be determined 
with the aid of the momentum theorem even for this case. This 
calculation was made by prandtl and Von Kcrmari.* The thickness 
of the boundary layer w a s found to be 

6 = 0.370 R x ~ l/5 (4) 



The shearing stress was 

U 6N 



o 



t 0 = 0.0225 p U 2 < 5 ) 



♦Zeitschrift fur angewandte Mathematik unci. Mechanik, Vol. I, 
1921, pp. 233-298. 
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The fundamental velocity distribution for 

/V\ l/7 

0<y<6 is u = U (jj- ; and for y > 6 is u =- U. 

^o j 

II. Researches of Burgers ?,nd Z.ijnan 

In 1923-1935, J. M. Burgers and B. G. van der Kegge Zijnen, 
in the Aerodynamical Institute of the Delft Technical High 
School, investigated the velocity distribution in the boundary 
layer for both kinds of flow and particularly in the region of 
transition from the laminar to the turbulent state.* In the 
laminar region they found very good agreement between the theo- 
retical and experimental results, as regards the thickness of 
the boundary layer. The discrepancies were greater, however, 
as regards the shearing stresses, whether determined from the 
velocity gradient on the plate, or from the momentum theorem 
with the help of the formula 

d 6 

T ° = d x P / U (U ~ u) d y (6) 

The . point of transition from the laminar to the turbulent region 
could be quite accurately determined by plotting the measured 
boundary-layer thickness as a function of the coordinate x. 
The resulting curve has a sharp bend at a certain value of R x . 
This value of R x may be designated as the "critical character- 
istic 11 (Konnzahl)o On the other hand, a "characteristic 11 can 

be aligned La the bevr $ a";y layer itsel f, b y express ing it jln 

* il Measurements of the Velocity Distribution in the Boundary Lay- 
er Along a Plane Surface," Delft, 1934. 
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the form Rg = Tj6/i?. Rg and R x are theoretically connected by 
the expression R5 = 5.5 J R x , which was also approximately 
confirmed by experiment. The measurements showed that R§ var- 
ies between 1650 and 3500" and is therefore considerably greater 
than in tubes. In the turbulent region the measurements showed 
that the velocity varied approximately as \fy. The curve of 
the boundary-layer thickness does not correspond to the formula 
6= 0,37 x R x ~ l/5 > because in its derivation, it was assumed 
that the turbulent flow begins at x = o, though, in reality, 
the laminar, flow continues up to the critical value of R x . 
Zijnen extended the formula by introducing a parameter x 0 and 

putting , 

6 - 0.37 (x - x 0 ) l/S (>) - (V) 

in which x 0 depends particularly on the shape of the leading 
edge of the plate ( shaxp or rounded). The experimental results 
can be represented very well by this modified formula. 

III. Object of Recent Experiments - Laboratory Equipment 

My own experiments were intended, on the one hand, to ac- 
count for the discrepancies in the above-mentioned experiments 
and, on the other hand," to determine the velocity distribution 
on rough plates. The means for these experiments were placed 
at my disposal by the "Notgemeinschaf t der Deutschen Wissen- 
schaft" (Emergency Fund for the Promotion of German Science), 

The experiments were performed in the small wind tunnel 
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belonging to the Aerodynamic Institute of the Aachen Technical 
High School (Figs, 1-2). The diameter of the entrance cone is 
30 cm (11.8 in.). The length of the free jet between the entrance 
and exit cones is about 2.5 m (8.2 ft.)* The measurements were 
made in this free jet where the static pressure was constant, 
which was essential for the method of measurement used. 

The wind velocity U was accurately controlled by enlarg- 
ing or diminishing the exit opening a by the use of plates* 
It was not feasible to regulate the wind velocity by means of 
a slide valve at a point b of the, enclosed portion of the wind 
tunnel, because this disturbed the air flow too much. This 
method was used, however, to produce artificial turbulence. 
Figure 3 shows the plates used in the tests, together with their 
degrees of roughness. The plates were successively mounted in 
a wooden frame on the test stand. Underneath the plate and 
mounted on a block of concrete, there was a lathe bed which sup- 
ported the measuring instrument. This arrangement made it pos- 
sible to move the measuring instrument either parallel or per- 
pendicular to the plate. x was read on a scale parallel to 
the plate. The perpendicular distance y from the surface of 
the plate was measured to within 0.01 mm (.0004 in.) by a slide 
gauge mounted on the same support. The pressure was determined 
by means of two alcohol pressure gauges. The temperature and 
\ pressure were read frequently during the tests. 

Though the Delft measurements were made with a hot-wire 
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instrument, I decided to use small Pitot tubes, which are com- 
paratively easy to make by drawing out and bending the ends of 
very small glass tubes. The smallest tube drawn by me had an 
outside diameter of 0.135 mm (0.0053 in.). They worked so slow- 
ly, however, that the readings would have taken too much time. 

Tube Sfo. 13, with an outside diameter of 0.35 mm (0.0138 
in.) and an inside diameter of 0c21 mm (0*0083 in.), was the 
most satisfactory and was used for most of the measurements. 
Figure 4 shows the calibration curves of this tube. The re- 
sults obtained with it were then compared with those obtained 
with an ordinary Pitot tube. If discrepancies were found, they 
could be mostly eliminated by grinding the open end of the tube. 
The measurements could be made only in a velocity region of 4 
to 36 m/s (13 to 118 ft. /sec), since the tubes vibrated 
strongly at higher speeds. 

In order to determine the distance between a tube and the 
plate, I" had to find the position at which the tube touched the 
plate. This position was easily determined by observing when 
the tip of the tube and its image came together (Fig. 5). 

A sufficiently accurate adjustment of the plate was effect- 
ed by mounting it perpendicular to the plane of the Pitot tube 
with the aid of a carpenter 1 s square. Any more accurate adjust- 
ment did not appear necessary, since measurements on both sides 
of the plate showed that the results were hardly affected by a 
•slight obliqueness of the plate. 
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IV. Results 



a) Laminar Flow 



Figure 6 shows the results for the thin dural plate No. 1. 
It is seen that the measured velocity distribution' agrees very 
well with the theoretically calculated velocity. The velocity 



ratio is plotted as a function of the quantity y/J and - 
shows that, within the accuracy of the measurements, the sec- 
tions very well obey the law of similarity. 

In Figure 7 the thickness of the boundary layer is plotted 



with the theory. 

For the thicker plate No. 4 (Figs. 8-9), the first sections, 
with small values of x, show systematic deviations which dis- 
appear for thicker sections. 

Figure 10 gives a verification of the impulse balance by 
comparing the shearing stress, as determined from the velocity 
gradient, with the value obtained from the momentum by means of 
formula (6). As shown in the figure, the two values obtained 
from my experiments agree very well, while there is a discrep- 
ancy in the values obtained from the measurements by Burgers. 
In my experiments there are differences at the leading edge of 
the plate. Hence it may be inferred that the slight discrepan- 
cies in my measurements are due to the influence of the leading 
edge, while there must be some other reason for the discrepant 




as a function of the quantity 




and likewise agrees well 
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cies in the measurements made by Burgers* 

In order to ascertain the effect of the finite thickness 
of the plate, especially the difference between the rounded and 
sharpened plates, I tried first to calculate the potential flow 
and then the curve of the boundary-layer thickness for such 
plates. The potential flow was calculated by the so-called in- 
verse method, by assuming source distributions which, when super- 
posed on the parallel flow, yielded streamline shapes similar 
to the plates used by me. The rounded section can be generated 
by pointlike sources, while distributed sources must always be 
adopted if the sharp section is to be substituted. 

Figure 11 shows the velocity curve calculated from the 
source distribution. In addition to the velocity curve, its 
first and second differential quotients are given according to 
the length of the arc, because these quantities are necessary 
to determine the thickness of the boundary layer. This thick- 
ness w a s calculated by A. Pohlhausen's method (Abh. Aerodyn. 
Inst. Aachen, No. l). Figure 12 shows the course of the quantity 
z = 6 3 /u for both the cases represented in Figure 11. With the 
rounded plate (pointlike source) the thickness of the boundary 
layer increases very rapidly and separation soon occurs. With 
the sharp-edged plate (linear source distribution), on the other 
hand, the deviation from the Blasius case (infinitely thin plate) 
is not very great. This result is in accord with the observed 
fact that, even at a short distance from the leading edge of the 
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rounded plate, the boundary layer becomes turbulent, while with 
the sharp-edged plate the laminar flow persists much longer. 

As regards the experiments of Burgers and Zijnen, this cal- 
culation could not explain the great discrepancy between the 
theoretical and experimental results as found, e.g., in the 
velocity gradient on the wall of the tunnel. The calculation 
might be utilized at most to explain the S-shaped velocity 
curve found by the above-mentioned investigators. The "boundary- 
layer theory, in fact, gives such curves, provided a pressure 
increase or velocity decrease (U 1 = -g— ^ < 0 in Fig. 11) is 
present. 

The reason why Burgers and Zijnen found a considerable de- 
viation from the theory of Blasius, while my measurements show 
a much better agreement, is due, in my opinion, to the' fact that 
the Delft experiments were performed in a closed wind tunnel, 
while mine were made in a free air stream. In a closed wind 
tunnel, a. decrease in pressure or an increase in velocity, cor- 
responding to. the frictional resistance of the tunnel walls 
and of the plate itself, must take place outside the boundary 
layer. The air flowing between the tunnel wall and the test 
plate is in a similar state to that in the first part of a tube. 
The resulting deviations from the Blasius theory can be approx- 
imately calculated, since. Zijnen himself established the increase 
in velocity along the test plate, though for another purpose 
(Thesis, Report Ho. 6, Delft, 1924, pp. 39-42). He expressed the 
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mean velocity increment by = f- U. The factor B denotes 

a x • l 

the velocity increment for the entire length I as a fraction 
of the velocity U. For example,' Zijnen found that 6 = about 
0.06 for U = 8 m/s (26 ft. /sec). I have investigated the 
effect of the velocity increment on the thickness of the bounda- 
ry layer according to K. Pohlhausen 1 s method and found smaller 
values for the thickness than the Blasius theory would lead one 
to expect. In fact, the Delft measurements show a deviation 
of the same nature. The velocity gradients on the plate can be 
compared for the two cases. For this purpose, in the calcula- 
tion of the boundary-layer thickness, I introduced ZijnerJs 
experimentally found values for 6 and the above value for 

U* = 4*— into Pohlhausen 1 s differential equation for the boun- 
d x 

dary layer. The result is shown in Figure 13. The plain lines 
represent the velocity gradient according to the Blasius method, 
while the dash lines represent the sa,me according to rny own 
method, with allowance for the velocity gradient in the tunnel. 
The agreement is very good. up to the last value, which corre- 
sponds to a velocity curve at the distance x = 62.5 cm (24.6 
in.) from the leading edge of the plate. With this velocity 
curve the transition to the turbulent condition has evidently 
taken place already, so that agreement can no longer be expected. 

The above-mentioned calculation also shows that the boundary- 
layer formation is very sensitive to relatively slight local 
variations in the velocity of the air flow along the plate so 
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that such measurements must really be linked with a very, accurate 
control of the static pressure. Thus, even an undulating or 
wavy shape of the plate may affect the formation of the laminar 
"boundary layer. 

b) ' Transition from Laminar to Turbulent Flow 

We have already remarked that the boundary layer can be 



assigned a critical characteristic value R§ = 5.5 J R x , at 
which the transition from the laminar to the turbulent condition 
takes place. I have attempted to determine R5 by two methods. 
The transition is evidenced, on the one hand, by a sudden growth 
of the boundary layer and, on the other hand, by the development 
of a shearing stress. Zijnen determined R5 by plotting the 
boundary- layer thickness. In a similar manner (Fig* 14) I first 
plotted the nondimensional quantity 6/ J and then also the 

quantity — , t— as a function of the characteristic 

U 3/2 (v x)" 1/2 

value. Both quantities must be constant for a purely laminar 
boundary layer. 

Despite the scattering of the test points, a bend is quite 
clearly indicated at R5 = about 3100. This result agrees with 
the Delft measurements. If this value of R5 is compared with 
the critical characteristic of tubes or rings, it appears very 
high. If it is considered, for example, that the thickness of 
the boundary layer should logically be taken as the hydraulic 
radius in our case, R5 would have a value about six times as 
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large as the corresponding value for tubes. It should be noted, 
however, that, even for tubes and rings, the value of the crit- 
ical characteristic can be greatly increased by a steady inflow. 
In the case of the sharp-edged plate, we have only a compara- 
tively steady flow. If the flow is disturbed, e.g., by the sep- 
aration of the boundary layer on a rounded leading edge, the 
transition takes place at a lower characteristic value. Also 
when the ear stream itself, for ex^jnple, due to strong throt- 
tling, was more turbulent than under normal conditions, I found 
that the reversal occurred sooner. The lower limit of R5 has 
not yet been established. Figures 15-16 give examples of ve- 
locity distribution in the transition region. 

o-l) Turbulent Flow over a Smooth Plate 

Since the preliminary experiments showed that the transi- 
tion to the turbulent flow is accelerated by rounding the lead- 
ing edge of the plate, a rounded strip of wood was applied to 
the front edge of the glass plate. With this arrangement a tur- 
bulent flow had already developed at x = 10 cm (3.94 in.). The 
tests covered a velocity range of 16 to 36 m/s (52-118 ft. /se.c.). 
Between x = 10 cm (3.94 in.) and 50 cm (19.69 in. measure- 
ments were made on the plate section at intervals of 5 cm (1.97 
in.). The corresponding shearing stresses on the wall were de- 
termined from formula (6). No direct determination of the shear- 
ing stresses on the wall could be made with the help of the ve- 
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locity gradient, because the measurement with the Pitot tube 
did not give sufficiently accurate results close to the wall. 
The results of the measurements are shown in Figure 17, both the 
distance from the wall and the velocity being plotted logarith- 
mically. It is obvious that, disregarding the immediate proxim- 
ity of the wall (y < 1 mm), the velocity distribution can be 
quite well approximated by a potential law. 

In Figure 18, the velocity ratio u/U is likewise loga- 
rithmically represented as a function of the ratio y/6 . The 
values of 6 were derived from the points of intersection of 
the sloping lines in Figure 17 with the horizontal lines in u/U. 
In other words, the lines corresponding to the different values 
of x were brought into coincidence by shifting horizontally. 
Figure 18 shows that the slope of the lines, i.e., the exponent 
of the potential expression depends somewhat on the velocity. 
If we put 

7? = o ( P (c = proportionality factor) (8) 

the exponent n then increases with the velocity. For example: 

n = 0.186 for U = 20 m (65.6 ft.); 
n = 0.196 » U = 28 « (91.9 « ); 
n = 0.198 " U - 36 " (118.0 « ). 

The exponents obtained for the smooth plate by this method are 
all higher than the exponent n = l/7 = 0.143- derived by Prandtl 
and von Karman from the law of resistance. This is due in .part 
to the fact that the test points y < 1 mm were disregarded in 
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in the evaluation. 

If the velocity distribution is approximated by a poten- 
tial law, the integral in formula (6) can then be evaluated. 
If the boundary-layer thickness is introduced in the above- 
mentioned sense, it is found that 

, 0 = fAapU* (9) 

in which the quantity a is connected with the exponent n by 
the expression 

. a = — I 1 , 

n + 1 2 n + 1 

I have usually determined the shearing stress with the aid of 
formula (9). 

If it is assumed that, between the shearing stress, the 
distance from the tunnel wall and the velocity, there exists a 
definite relation governed only by the physical constants, den- 
sity and viscosity, this relation can then be expressed in the 
form 

= f f^y\ do) 

p U S V V / 

This equation is the usual mathematical expression of the Prandtl 
dimensional consideration, which leads to the previously men- 
tioned formula for the turbulent velocity distribution. The 
quantity uy/u can be regarded as a kind of characteristic coef- 
ficient R y for the point y.. If the velocity distribution is 
approximated by a potential formula, equation (10) is changed to 
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T in El 

p u 2 \ v J 

The relation between the exponents n and m is evidently ex- 

pressed by the formula m = - n"T~T # ^ is a coe:f f icien *k 
which can "be determined from the measurements, 

Theodor von KaJfman, on the basis of the Blasius law of re- 
sistance, arrived at the 'values: t = 0.0235, m = 0.25, 
n = 0.143. I have compared all the data obtained for the smooth 
plate with the von Karman theory (Fig. 19). It is seen that 
the points for small values of uy/v , . in the vicinity of 
u y/v g.2000, scatter about the theoretical curve, which a sys- 
tematic deviation occurs for greater values. Greater values 
than n = 0.143 are likewise obtained by disregarding the val- 
ues in the immediate vicinity of the tunnel wall, as already 
mentioned. 

c-2) Turbulent Flow along Rough Plates 

After the measurements on the smooth plate, I determined 
the velocity distribution on an undulating plate and on two 
rough plates. It was found that the interpolation with the po- 
tential formula gave good results even in this case. The data 
were worked out in the same manner as for the smooth plate 
(Figs. 20-22). 

Furthermore, the resistance and the shearing stress were 
determined from the impulse integral. For comparison I deter- 
mined the integral first on the basis of the test points, but 
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then also with the aid of formula (9). The deviations were 
greater for the rougher plate, apparently because the interpola- 
tion with the aid of the potential law does not hold good up. 
to values of u = U or y = 6, but the velocity-distribution 
curve is sharper. In this case the exact integral was naturally 
resorted to for the calculation of the shearing stress. 

The next step was to determine to what extent the shearing 
stress can be represented by formula (11). Figures 23 and 24 
show the result. The quantity T 0 /p u a is plotted as the ordi- 
nate and uy/v as the abscissa, both in the logarithmic scale. 
It is obvious that the absolute value of the exponent m in- 
creases, on the one hand, with increasing roughness and, on 
the other hand, with increasing velocity. 

For all the data, the coefficient £ varies between 0.03 
and 0.6, the exponent m between -0.3 and -0.5 and correspond- 
ingly the exponent n between 0.186 and 0.325. I have now made 
the noteworthy observation that the values of £, as a function 
of n in the logarithmic scale, lie on a straight line. 

In my opinion, the fundamental importance of the results 
resides in the fact that there is apparently a definite rela- 
tion- between the exponent of the velocity distribution and the 
coefficient £ . The latter corresponds to the coefficient which 
is assumed to be constant in the theory for the ideal case of 
the smooth plate . For rough plates, however, both the quantity 
£ and the exponent m are no longer constants, but are affected 
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by the velocity as well as by the roughness. It may be surmised 
that both depend on the nondirnensional quantity uk/u , in which 
k is the mean roughness. The general resistance formula 
would then have the form 

Jlo^ = C±J£\ K V 7 > f3U (13) 
p u 3 V v / ' V v / 

It is still uncertain whether this equation represents only an 
interpolation formula or a definite law. 

' Translation by 
National Advisory Committee 
for Aeronautics. 
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Figs. 11,12 
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Figs. 13,14 




_„ // Burgers % / / Burgers 
S // , p §/ 

o */7 




Blasius ^ / Blasius- 



y x=10 cm 



x=20 cm 



Burgers 



Burgers 




Burgers 



x=30 cm x=40 cm x=62.b cm 
Fie. 13 



6^ 
vx 
U 



o 



CV2 



-.16 
14 
12 
10 
8 
6 



i to 4 
^ I 
X o 

P H 

2 

— -H 

to 
to 





\ 














! 

























/ 






- 1 
_J 














— T^7 














„. 




r i/i 


















1 


'A" 






















+ 








































+ 


/ 










































+ 
























A- 




















?A 






















1 

'+ + 












































































4 






M' 
























4 ■ l— 


14- 


! 




H 

I 





From boundary- 
layer thick- 
ness 6 



From shearing 
stress Tq 

a, 1000] 

b, ?22^ v alue8 of R6 

c, 3000 

d, 4000 



a Id c" R x in 10^ 



Fig. 14 



N . A . C . A . T e clrni cal Memorandum No . 585 



Figs. 15,16 



1 


* 






1 

X 

\ 

1 1 
I 


c 

\ 

X 


Plato 4 


= 16 m/s 


X 

X 

\ 




*-x- 




1 




"*-x-^ 




X 








X 


*--x- 


— x 




i 

X 








it 




-x 




t 


Mil 


ill! 


Mil 



LO 



O LO 



o 

O' 

LO 



o 
to 



o 

H 

LO> 



! 








> 


< 






CQ 




1 






a 


i 

> 
> 


X 

c x 
l x 

X 

X X 

5 N 

** 

X 


*x 


o 

1 -> 

H 
Ph 

' 


U = 20 




X 

X 




* 






< 

X 












*' 


-vc 






111! 


Mir 


Sill 


MM 



c A 



a 

CO* 



02 



O 
CQ 



LO 
H 



O 
H 



LO 



o 

LO 

.ft 



H 
II 

a 

o 

H 



LO 
H 

bb 

•H 



S/tu c Vi 



\ 

X 

\ 

x x 

\ \ 
\ 



W y 

II I 1 I X 



to 
-p 

H 



*~X-y 



I ( I X 



CQ 

a 

CD 

H 
II 



I 1 IX 



A 



a 

O 
LO 

o 

LO 



o 

H 



a 

o 



LO 
H 



O 



LO 



s /iu f n 



H 

LO >> 

CQ 



\ 

X 

\ 

\ 

X X 
X X 
X X 

w 

X X 



JJJ_L 



IXXX 





CQ 






CD 




■P 


cv? 


cd 




H 


II 


PL, 









MM 



Mil 



LO 
CQ 



O 

cq 



LO O LO 

H rH 

s /m ' n 



CO 
H 

•H 



a 
o 

LO 



LO 



a 

b 



o 

LO 



